The Effects of Asymmetric Salt and a Cylindrical Macroion on Charge Inversion: 
Electrophoresis by Molecular Dynamics Simulations 
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The charge inversion phenomenon is studied by molecular dynamics simulations, focusing on size 
and valence asymmetric salts, and a threshold of surface charge density for charge inversion. The 
charge inversion criteria by the electrophoretic mobility and the radial distribution functions of 
ions coincide except around the charge inversion threshold. The reversed electrophoretic mobility 
increases with the ratio of coion to counterion radii a" /a + , while it decreases with the ratio of coion 
to counterion valences Z~ /Z + . The monovalent salt enhances charge inversion of a strongly charged 
macroion at small salt ionic strength, while it reduces reversed mobility otherwise. A cylindrical 
macroion is more persistent to monovalent salt than a spherical macroion of the same radius and 
surface charge density. 

PACS numbers: 61.25.Hq, 82.45.-h, 82.20.Wt 



I. INTRODUCTION 

The phenomenon of reverting the charge sign of large 
ions due to other ions and salts in water solution was 
known to physical chemists as charge inversion or over- 
screening for half a century More recently, it also 
attracted a si gnif icant attention of physicists 0, H H IE 
SUSS MM El El El El El Gj Charge inver- 
sion is now understood to be the generic phenomenon 
that occurs in strongly correlated charged systems. It 
has far reaching consequences in biological and chemical 
worlds. In particular, it seems to be a decisive ingredient 
in modern gene therapy, facilitating the delivery of genes 
(negatively charged DNA) through cell walls of predom- 
inantly negative potentials 0, El . 

In our previous papers, we adopted the static E3 an d 
dynamical |13| models to study charge inversion of the 
spherical macroion by molecular dynamics simulations. 
In the former study where a macroion was immovable, 
the static quantities were obtained including the radial 
distribution profiles of counterions and coions moving in 
the Langevin thermostat. Generally speaking, the ra- 
dial profiles of charge density and integrated charge are 
good indices of charge inversion. However, it is not so 
when the charge inversion is around the threshold (a 
counterexample is Fig^b) of this paper). Namely, those 
ions forming the outskirt of the radial distribution profile 
are not electrostatically bound to the macroion complex 
and are left behind when the complex moves in the so- 
lution. To separate the bound ions from the floating 
ones, the electrophoretic mobility problem in the explicit 
(particle) solvent was studied in the latter. We showed 
that a macroion satisfying the charge inversion conditions 
drifted with adsorbed counterions and coions toward the 
external electric field, the direction of which agreed with 
the sign of inverted charge. 

This paper adopts the dynamical (electrophoresis) 
model and focus on the following points. First, we ad- 
dress the asymmetry effects of the size and valence of 
counterions and coions. The radii of bare halogen atoms 



are large compared to those of alkaline metals, which be- 
come eventually larger than the former after hydration. 
Secondly, we examine the charge inversion threshold in 
terms of the surface charge density and monovalent salt 
for the macroions of spherical and cylindrical shapes |20| . 
The adsorption of cylindrical macroions to a charged sur- 
face was previously studied H 0] . Instead, we deal with 
the charge inversion of a cylindrical macroion by adsorp- 
tion of counterions in the presence of monovalent salt. 
The cylindrical macroion is shown to be more persistent 
to monovalent salt than the spherical macroion of the 
same radius and surface c harg e density. 

In the dynamical study [l3( , the average static quanti- 
ties including the radial distribution profiles of ions and 
the structure of the electrical double layer are not altered 
in the zeroth order as far as the applied electric field E 
is much smaller than the electric field produced by the 
macroion charge, E <C Qo/eRg, which is a huge electric 
field 10 6 V/cm realized by high-power short pulse lasers. 
Here, Qo and Rq are macroion bare charge and radius, 
respectively. We have confirmed this by excellent agree- 
ment of the integrated charge profiles between the runs 
with and without the applied electric field, in 3% accu- 
racy for equal-sized coions and trivalent counterions with 
the parameters of Fig|5]in Sec. 3. The net charge of the 
macroion complex Q* can be estimated if one makes use 
of the force balance relation, 



Q*E-vV^0, or Q* 



(1) 



where \l = V/E is the measured electrophoretic mobil- 
ity and v is the solvent friction. However, one needs 
to assume the friction for this conversion. The fric- 
tion that acts on the charged object in charge-neutral 
electrolyte solution is an enhanced one compared to the 
Stokesian friction because the velocity decay length is 
drastically shortened by screening of hydrodynamic inter- 
actions |2lj, |22( . Despite of this uncertainty, the observed 
mobility agreed qualitatively well with the net charge 
scaling derived by theory pJJ | . 

Other issue of the dynamical study is that the Joule 
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heat is produced by the applied electric field, which is 
transferred to background neutral particles through col- 
lisions with the ions. To drain the heat, we adopt a heat 
bath at the domain boundary that is kept at a constant 
temperature T. The thermal bath screens hydrodynamic 
interactions and nullifies the momentum of the system on 
time average. However, it is emphasized that the electric 
field does not bring in momentum into the charge-neutral 
system, and that hydrodynamic interactions are weak 
compared to the electrostatic binding forces and screened 
at short distances in the electrolyte solution [2ll . l22T | . Mo- 
mentum is not transported beyond this distance from the 
macroion toward the domain boundaries. We confirmed 
this fact by comparing the mobilities of two runs with 
and without the heat bath (before heating became sig- 
nificant), the results of whic h ag reed very well within 
simulation errors (cf. Fig. 5 of [l3). Thus, the heat bath 
does not affect our simulations. 

This paper is organized as follows. The simulation 
method and parameters are described in Sec[Hl The ef- 
fects of asymmetric radii and valences of coions and coun- 
terions on the charge inversion phenomenon are stud- 
ied in Sec lIIII The dynamical and static observables, 
i.e. electrophoretic mobility and radial distribution func- 
tions of ions, are compared for the same runs, where 
good correspondences are found except around the charge 
inversion threshold. In Sec II VI the effects of monova- 
lent salt that exits as the base of the Z:l counterions 
and coions are examined. There, a cylindrical infinite 
macroion (occupying full length across the domain) that 
can adsorb geometrically more counterions than a spheri- 
cal macroion of the same radius is introduced. The cylin- 
drical macroion is found to be more persistent to mono- 
valent salt than the spherical macroion of the same radius 
and surface charge density. SeclVI will be a summary of 
this paper. 

II. SIMULATION METHOD AND 
PARAMETERS 

Simulation method and parameters for the present 
study are described here. We take the system consist- 
ing of one macroion, many counterions, coions and neu- 
tral particles. We solve the Newton equations of motion 
for each particle with the Coulombic and Lennard- Jones 
potential forces under a uniform applied electric field E 
(E > 0). A large number of neutral particles are used to 
model the solvent of given temperature and to treat the 
interactions among the finite-size macroion, counterions 
and coions. 

The units of length, charge and mass are, a, e, and to, 
respectively. As is mentioned below, our choice of the 
temperature corresponds to a ~ 1.4A in water and to ~ 
40 a.m.u. A spherical macroion with radius Rq = 5a, 
negative charge Qq between — 15e and — 81e, and mass 
200to is surrounded by the N + number of counterions of 
a positive charge Z + e and the N~ coions of a negative 
charge —Z~e. The corresponding surface charge density 



of the macroion is between a sp = 0.048e/a 2 (0.39 C/m 2 ) 
and 0.26e/a 2 (2.1 C/m 2 ). A rod-shaped macroion is also 
used in Sec lIVI whose surface charge density is between 
(T rod = 0.04e/a 2 (0.33 C/m 2 ) and 0.08e/a 2 (0.66 C/m 2 ). 
The system is maintained in overall charge neutrality, 
Qo + N + Z + e — N~ Z~e = 0. The radii of counterions and 
coions are a + and a~, respectively, with the counterion 
radius being fixed at a + = a, and the radius of neutral 
particles is a/2. The mass of the coions and counterions 
is to, and that of A* neutral particles is to/2, where we 
have in mind the water molecule against that of K + or 
Ca 2+ ions. However, the mass is not involved as far as 
an equilibrium state is concerned. Approximately one 
neutral particle is distributed in every volume element 
(2.1a) 3 w (3A) 3 inside the simulation domain, excluding 
the locations already occupied by ions, which typically 
yields 4000 neutral particles for the spherical macroion 
of radius i?o = 5a. These particles are placed in a cubic 
box of size L = 32a, with periodic boundary conditions 
in all three directions. 

Calculation of the Coulomb forces under the periodic 
boundary conditions requires the charge sum in the first 
Brillouin zone and their infinite number of mirror images 
(the Ewald sum 23]). The sum is calculated with the use 
of the PPPM algorithm [25j|. We use (32) 3 spatial 
meshes for the calculation of the reciprocal space con- 
tributions to the Coulomb forces, with the choice of the 
Ewald parameter a « 0.262 and the cutoff length 10a. 

In addition to the Coulomb forces, the volume exclu- 
sion effect among particles is introduced through the re- 
pulsive Lennard- Jones potential (/>lj = 4e[(A/rij) 12 — 
[A/r l3 f] for m = |r ( - tA < 2^ 6 A, and </> LJ = -e 
otherwise to exclude the attraction part. Here is the 
position vector of the z-th particle, and A is the sum of 
the radii of two interacting particles. We relate e with 
the temperature by e = fc^T, and choose ksT = e 2 /5ea 
(we assume spatially homogeneous dielectric constant e). 
The Bjerrum length is thus = e 2 /efcgT = 5a, which 
implies a w 1.4A in water. 

The initial states of the runs are prepared by randomly 
positioning all the ions and neutral particles in the simu- 
lation domain and giving them Maxwell-distributed ran- 
dom velocities. The Newton equations of motion are in- 
tegrated with the use of the leapfrog method, which is 
equivalent to the Verlet algorithm |26J . The unit of time 
is r = ay/m/e, and we choose the integration time step 
At = 0.0 It, where t s» lps in real environments. The 
simulation runs are executed well beyond the time when 
the dynamical property (the macroion drift speed in this 
study) has become stationary, and only the stationary 
part is used for data analysis. Under the electrostatic 
and Lennard-Jones forces where the random forces and 
frictional forces are represented by particle solvent, typ- 
ical run times are between 2000t and 4000r (some runs 
are continued further to check the long term variability, 
but the results fall within the error bars). 

As noted earlier, we can adopt the heat bath with- 
out side effects for the present study with the neu- 
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tral electrolyte solvent. At the center of the heat bath 
the macroion is located at every moment. The veloci- 
ties of the neutral particles crossing the boundaries are 
refreshed according to the thermal distribution of the 
given temperature. Throughout this paper, the same 
normalization of the mobility n$ = «o/(|<3oo|/e-Ro) ~ 
21(/zm/sec)/(V/cm) is used, where Q a — — 80e, R = 5a 
and vq is the thermal speed of neutral (solvent) particles. 

III. THE EFFECTS OF ASYMMETRIC 
COIONS/COUNTERIONS 

The dependence of the electrophoretic mobility of a 
macroion on the ratio of coion and counterion radii 
a~/a + is shown in Fig^ The coion radius is varied 
between 0.5a and 2.5a for the fixed counterion radius 
a + = a, and the Bjerrum length is As = e 2 /eksT = 
5a. The counterion valence is Z + = 2 (diamonds), 3 
(circles) and 5 (squares). The number of monovalent 
coions is N~ ~ 60, which corresponds to the density 
1.8 x 10~ 3 /a 3 wl.l Mol/L and that of the counterions 
is determined by charge neutrality. The macroion with 
charge Qo = — 80e and radius Rq = 5a has the surface 
charge density a sp ~ 0.26e/a 2 (2.1 C/m 2 ). The electric 
field is E = O.le/ae -C Qo/eR$ which ensures that the 
electrophoresis occurs in the linear regime [Tflj . 

Figure ^ shows that the mobility is reversed (posi- 
tive) and increases almost linearly with the ratio a"/a + . 
The mobility peaks at the intermediate coion radius 
a~/a + 1.5, beyond which the mobility sharply de- 
creases and flips to non-reversed. The degradation of the 
mobility does not depend on the counterion valence nor 
is indexed to the coion and counterion interaction energy 
Z + Z~e 2 /e(a~ + a + ) which is always larger than ther- 
mal energy ksT (between 5 and 10 times for Z + = 3). 
On the other hand, the volume fraction of the ions be- 
comes significant, 47r(7V-(a-) 3 + N+ (a+) 3 )/3L 3 - 0.07 
for a~ = 2a and 0.13 for a~ = 2.5a. 

The reversed mobility at a~ = a + and Z + = 3 is 12 
(/im/sec)/(V/cm). The experiment with a sulphonated 
latex particle of the surface charge density 0.115 C/m 2 
obtained the saturated mobility 0.25 (um/sec)/(V/cm) 
for 50 mM/1 La(N0 3 ) 3 (Fig. 6 of Ref.^17]). Although ex- 
act comparison is difficult, we note that the theory pre- 
dicts [l(| that net charge Q* , hence the mobility, scales 
as Q* /R 2 , ~ a(Rws I '-Vd) ~ (Zeacg) 1 / 2 for the weak De- 
bye shielding regime, where A_d = (efc^T '/Snoge 2 ) 1 / 2 with 
c s the coion density, and Rws = (Ze/nasp) 1 ^ 2 is the 
Wigner-Seitz cell radius. This formula may be applica- 
ble because of weak lateral screening on the macroion 
surface due to the absence of coions there. The mobility 
of the former is expected to be 37 times that of the latter, 
whose agreement is fair. 

The cases with reversed mobility are characterized by 
the peaked radial distribution profile of the counterions 
localized on the macroion surface and strong association 
of coions to the counterions, as seen in FigE^a). The 
integrated charge Q(r) in the inset panel is peaked and 
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FIG. 1: The electrophoretic mobility of the macroion /i is 
shown against the ratio of coion and counterion radii a~ /a + . 
The counterion radius is fixed at a + — a, and their valence 
is Z + — 2 (diamonds), 3 (circles), and 5 (squares). The sur- 
face charge density of the macroion is cr sp ~ 0.26e/a 2 (2.1 
C/m 2 ) for the charge Qo = — 80e and radius Ro = 5a. Here, 
the normalization of the mobility is fio = «o/(|<3o|A-Ro) ~ 
21(/^m/sec)/(V/cm) with Vo the thermal speed of neutral 
particles. The number of monovalent coions is TV - ~ 60 
(l.lMol/1). The Bjerrum length is A s = 5a. 
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FIG. 2: The radial distribution functions of charge density p s 
of counterions (solid bars) and coions (shaded bars) , with the 
coion radius (a) a~ — a and (b) a~ = 2a, are shown for the 
runs in Fig0(the coion bars are enlarged by ten times). The 
counterion radius is a + — a and valence Z + — 3. The inset 
panels show the integrated charge distribution Q(r)/\Qo\ of 
the corresponding counterion and coion charge densities. 



larger than the neutrality |Qo|- By contrast, for the case 
of non-reversed mobility at a~ = 2a in Fig[2[b), the ag- 
gregates made of one counterion and a few condensed 
large coions distribute more homogeneously and apart 
from the macroion than for the a~ /a + — 1 case. The 
radial profile of the integrated charge is nearly flat due 
to less concentrated counterions, thus the dynamical and 
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static observables of charge inversion are consistent. 

The increase in the mobility with coion radius arises 
from less amount of charge compensation by large coions 
due to geometrical avoidance on the counterion surface. 
This is consistent with the Monte Carlo simulation of 
charge inversion for the finite-size coions Q and the con- 
densation of the Z : 1 ions with the size asymmetry |27| . 
Our observation shows that the degradation of the mobil- 
ity for large coions is due to collisions; a part of adsorbed 
counterions are kicked out of the macroion surface whose 
return is hampered and delayed by large coions surround- 
ing the macroion. This is contrary to the charge inversion 
enhancement by the excluded volume ordering by large 
and symmetric coions and counterions |28| . 

The effect of asymmetric valences between counteri- 
ons and coions on the electrophoretic mobility is shown 
in Fig|3 The coion valence Z~ is varied for the fixed 
counterion valence Z + which is either divalent, trivalent 
or tetravalent. The macroion surface charge density is 
a sp ~ 0.26e/a 2 (2.1 C/m 2 ) for the charge and radius 
Qo ~ — 80e and Ro — 5a. The charge content carried 
by coions is kept the same such that N~ = 300/Z", and 
that of the counterions is determined by charge neutral- 
ity. The radius of the coions and counterions is equal, 
a~ = a + — a. 

Figure |3[a) shows that the mobility is reversed and 
largest for the monovalent coions, i.e. for the largest 
asymmetry with any Z + value. The mobility becomes 
larger as the valence of the counterions increases from di- 
valent, trivalent to tetravalent at Z~ = 1. As the degree 
of the asymmetry decreases with the increase in coion 
valence, the magnitude of reversed mobility decreases as 
until the two valences become nearly symmetric, 
Z + ~ Z~ . The mobility remains reversed for the equal 
valence cases except for Z + = Z~ = 1, and finally be- 
comes non-reversed for larger coion valences. If FiglSJa) 
is replotted for the coion and counterion interaction en- 
ergy modified by a factor 7 as in Fig|3Jb), three curves of 
different valences Z + follow a master curve, where 7 = 
1, 2 and 2.5 for Z+ = 2, 3 and 4, respectively. The 7 
factor stands for the reduction of the interaction energy 
by enhanced screening possibly due to condensed coions 
to the counterions, which is actually observed. 

The charge inversion for divalent counterions and 
coions is in line with the theory result of the HNC-MSA 
integral equations |14| . The radial distribution profile of 
the integrated charge is peaked Q pe ak/Qa ~ 1-4, and its 
trailing tail oscillates periodically around the neutrality 
as found previously [12J, |28|. Although the Z + = Z~ 
case resembles the Debye screening in the sense that the 
ratio of the macroion, counterion and coion valences is 
Qo/Z : 1 : 1, neither the Debye or nonlinear Poisson- 
Boltzmann theory is applicable due to strong electro- 
static interactions and charge inversion takes place. 

The dynamical results of the mobility in Fig|3|are com- 
pared with the static results in Fig0|that depicts the ra- 
dial distribution functions for trivalent counterions and 
the coions with (a) Z~ = 1 and (b) Z~ = 3. In the 





FIG. 3: The macroion mobility is shown against (a) the coion 
valence Z~ , and (b) the modified interaction energy of coun- 
terions and coions Z + Z~ e 2 /2eja, for the fixed counterion va- 
lences, Z + = 2 (triangles), Z + — 3 (open circles), and Z + — 4 
(solid circles). The 7 factor is 1.0, 2.0 and 2.5 for Z + = 2, 3 
and 4, respectively. Coions and counterions are of the same 
radius, a~ — a + = a, and the temperature is e 2 /eafcsT = 5 
(the Bjerrum length As = 5a). 



panel (a), the counterions are strongly adsorbed to the 
macroion and the radial profile of the integrated charge 
in the inset panel is peaked, \Q pea k/Qo\ ~ 1-6. Both the 
dynamical and static observables indicate strong charge 
inversion. It is counterintuitive, however, that the max- 
imum peak height of the integrated charge in Fig0Jb) 
is not proportional to the small mobility measured in 
FigEl The integrated charge profile with a sharp dip 
after the first peak is usually associated with small re- 
versed or non-reversed mobilities. For crudely detecting 
the charge inversion statically, one has to count the num- 
ber of the ions whose binding potential to the complex is 
larger than fcgT, which is not trivial in data analysis of 
the molecular dynamics simulations. 



IV. SPHERICAL AND CYLINDRICAL 
MACROIONS UNDER MONOVALENT SALT 

We study in this section both the spherical and cylin- 
drical macroions under the monovalent salt that exists as 
the base component to the Z:l multivalent salt. Here, the 
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FIG. 4: The radial distribution functions of charge density p s 
of counterions (solid bars) and coions (shaded bars) , with the 
coion valence (a) Z~ — 1 and (b) Z~ = 3, are shown for the 
runs in Fig[3] The counterion valence is Z + = 3. The inset 
panels show the integrated charge distribution Q(r)/\Qo\ of 
the corresponding counterion and coion charge densities. 
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FIG. 5: The electrophoretic mobility of the spherical 
macroion is shown against ionic strength of monovalent salt 
n s i (1 Mol/1 salt corresponds to 0.0017/a 3 ). The surface 
charge density of the macroion is a sp — Qq/^tyR^ ~ 0.26e/a 2 , 
with excess Z-ions (filled squares), and without excess Z- 
ions (open squares). Also, the mobility for the macroion 
a sp ~ 0.080e/a 2 with excess Z-ions (filled circles) is shown. 
Here, fio ~ 21(/im/sec)/(V/cm). 



trivalent counterions (Z = 3) are called the Z-ions, and 
other counterions and coions are monovalent. All these 
ions have equal radius a~ = a + = a. 

Figure [3] shows the dependence of the mobility of a 
spherical macroion against ionic strength of monovalent 
salt, n s i = 2N +1 /L 3 , where N +1 is the number of mono- 
valent counterions (the same number of matching mono- 
valent coions are present). The N +1 = 50 salt ions cor- 
respond to n sI ~ 0.0031/a 3 (1.8 Mol/1), and the Debye 
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FIG. 6: The radial distribution functions of charge density 
p s of counterions (solid bars) and coions (shaded bars) are 
shown for (a) the spherical macroion in Fig|3 and (b) the 
cylindrical macroion in Fig|H] both with surface charge density 
a — 0.08e/a 2 , radius R — 5a, and salt ionic strength n s i ~ 
0.006/a 3 (the coion bars are enlarged by three times). The 
inset panels show the integrated charge distribution Q(r)/|Qo| 
of the corresponding counterion and coion charge densities. 



screening length is 1.6a. The surface charge density for 
the strongly charged macroion with Qo = — 81e (filled 
and open squares) is a sp ~ 0.26e/a 2 , and that for the 
weakly charged macroion with Qq = — 25e (circles) is 
0.08e/a 2 (0.65 C/m 2 ), both with the radius R = 5a. 
For these runs, the surface charge is at least several 
times that of the salt ions contained in the surface layer, 
a sp 3> en s i\]j. The Guoy-Chapman length is even 



smaller, Xqc 



counterions and a sp ~ 0.26e/a 2 . The Wigner-Seitz cell 



eksT '/2nZea S p 
T sp - 0.26e/a 2 . 
radius for the spherical macroion is, 

(Ze/Tra sp )^ 2 = 2R Q (Ze/\Q Q \) 1 / 2 



0.13a for trivalent 



R 



ws 



(2) 



This yields Rws ~ l-9a and 3.5a, respectively, for the 
strongly and weakly charged macroions in Fig[SJ A se- 
ries of the runs with spherical macroions of different bare 
charges Qo reveal that the charge inversion threshold un- 
der the zero monovalent salt is Qo ~ 15e for Rq = 5a. 
This yields Qo/R 2 , ~ 0.6e/a 2 and the surface charge den- 
sity at the threshold 



' sp 



0.048e/a 2 (0.39C/m 2 ). 



(3) 



This agrees with our previous finding [l3| that the thresh- 
old of surface charge density is \Qo\/Rq « 0.5e/a 2 for any 
combination of Qq and Rq . The correlation energy of the 
surface Z-ions forming the Wigner-Seitz cell becomes at 
the threshold 



Z 2 e 2 /2eR ws ~ 5fc B T. 



(4) 
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FIG. 7: The bird's-eye views of all ion species are shown for 
the spherical macroion with (top) and without (bottom) the 
excess Z-ions for the macroion of surface charge density a sp ~ 
0.26e/a 2 in Fig|S] The trivalent and monovalent counterions 
are drawn in red and yellow, respectively, and the monovalent 
coions in green (nearly 4000 solvent particles are not drawn). 



In FigO we see that the addition of small amount of 
monovalent salt enhances the reversed mobility for the 
strongly charged macroion with excess Z-ions N^~ 3 = 10 
(filled squares). The mobility peaks at n s j ~ 0.0015e/a 2 
(0.83 Mol/1). However, further addition of monovalent 
salt whose charge concentration exceeds that of the Z- 
ions in the bulk screens the electric field and suppresses 
the charge inversion. By extrapolation, the reversed 
mobility is expected to terminate at n s j ~ 0.013/a 3 
(7.2 Mol/1). Even without the excess Z-ions for which 
N +3 = \Qo\/eZ, the addition of monovalent salt induces 
the charge inversion (open squares; the fifth data point 
overlaps on the filled square). On the other hand, for 
a weakly charged macroion with excess trivalent Z-ions 
iV^ 3 = 10, the mobility is reversed at zero salt but de- 
creases monotonically with the salt ionic strength. 

A close look at the ion distributions shows that the en- 
hancement of the mobility at the small salt ionic strength 
in Fig0 is due to the rearrangement of ions to achieve 
global energy minimization rather than the adjustment 
of the Z-ions on the macroion surface. This is found in the 
radial distribution functions of ions in FigEfa), and visu- 
ally in the bird's-eye view plot in FigH For the case with 
excess Z-ions, almost all the Z-ions are adsorbed to the 



macroion, while monovalent salt ions including counteri- 
ons are located apart from the macroion surface. For the 
case without the excess Z-ions, both trivalent and mono- 
valent counterions are adsorbed to the macroion, but the 
coions are absent in the vicinity of the macroion. Thus, 
the lateral screening on the macroion surface should not 
be significant. The theory of charge inversion flfl ] pre- 
dicts enhancement of overcharging by monovalent salt. 
However, the theory heavily relies on the screening of 
surface Z-ion correlations beyond the nearest neighbor 
cells, and the predicted overcharging occurs in the lim- 
ited regime. 

We interpret that the aforementioned rapid decrease 
in the mobility with salt ionic strength for the weakly 
charged macroion is due to the finite radius of the spher- 
ical macroion. Namely, a number of Z-ions adsorbed on 
the small surface of the macroion may not be sufficient to 
maintain electrostatic correlations and overcome thermal 
fluctuations. 

Therefore, the rest of this paper is devoted to the 
charge inversion of a cylindrical macroion. Here, the 
macroion is assumed to be a rod of finite radius R ro d 
and infinite length, occupying the full length across the 
domain, and that it lies perpendicularly to the applied 
electric field for the ease of the simulation. It is remarked 
that a polyelectrolyte chain tends to align along the elec- 
tric field due to the polarization effect (2i|. However, 
the orientation of the rod may be of the second impor- 
tance in the case of a non-deformable rod macroion since 
the counterion adsorption by the macroion is not much 
altered as far as the applied electric field is weak com- 
pared to the local electric field produced by the macroion, 

rod^- 

Before other results are shown, the dependence of the 
mobility on the length of the rod is examined. Table I 
summarizes the mobility for the finite length rod which 
has the cylinder length i and is capped by two semi- 
spheres of radius R at each end. The surface area of 
the macroion is S = 2ttR(2R + £), and the fixed surface 
charge density o~ ro d ~ 0.08e/a 2 determines the macroion 
charge Q ro d- The I = case is a spherical macroion. The 
surface electric fields at the center and tips of the rod 
are thus almost the same for different rods. The mass 
of each rod is varied such that the ratio Q r od/M roc i w 
const, although the rod mass should not be involved in 
the equilibrium drift speed and the mobility. Table I 
shows that the mobility slightly increases for a short rod 
compared to the spherical macroion but becomes nearly 
insensitive to the length of the rod for £/R ro d > 3, if the 
surface charge density is maintained and the monovalent 
salt is not present. This rod length may indicate the 
transition from the 3-D to 2-D regime of the cylindrical 
macroion. 

The charges of the infinite rod macroions in Fig|S] are 
Qrod — — 80e (filled circles), — 60e (open circles) and 
— 40e (filled squares). The surface charge densities of 
these macroions are a ro d — \Qrod\/2TrR ro dL ~ 0.08e/a 2 , 
0.06e/a 2 and 0.04e/a 2 , respectively, which correspond to 
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TABLE I: The dependence of mobility fj, on the length £ of 
the cylindrical macroion with the radius R — 5a that lies per- 
pendicularly to the applied electric field. The surface charge 
density a roc i ~ 0.08e/a 2 gives the macroion charge Q ro d = cS 
for the surface area S — 2irR(2R + £). The number of coions 
is iV~ ~ 90 and monovalent salt is not present. 



i/a 





5 


10 


15 20 


S/a 2 


314 


471 


628 


785 942 


\Qrod\/e 


25 


38 


50 


63 75 




0.29 


0.35 


0.34 


0.31 0.30 



3L 




0.000 



0.002 



0.004 



0.006 



0.008 



FIG. 8: The electrophoretic mobility of the rod-shaped 
macroion is shown against ionic strength of monovalent salt 
n s j. The macroion is an infinite rod with the radius R ro d = 
5a. The surface charge density of the macroion is a ro d ~ 
0.08e/a 2 (filled circles), 0.06e/a 2 (open circles), and 0.04e/a 2 
(filled squares), which correspond to 0.66C/m 2 , 0.49C/m 2 
and 0.33C/m 2 , respectively. The normalization of the mo- 
bility is /io ~ 21(/xm/sec)/(V/cm). 




FIG. 9: The bird's-eye view of the cylindrical macroion and 
other ions located within 4a from the macroion surface are 
shown for the run in Fig|S| The surface charge density is 
a = 0.08e/a 2 (0.66 C/m 2 ) and the salt ionic strength is 
n s i ~ 0.006/a (3.6 Mol/1). The trivalent and monovalent 
counterions are drawn in red and yellow, respectively, and the 
monovalent coions in green (solvent particles are not drawn). 



0.66 C/m 2 , 0.49 C/m 2 and 0.33 C/m 2 . For reference, the 
surface charge density of the DNA is 0.19 C/m 2 . The 
number of Z-ions is adjusted as N +3 « \Q ro d\ /eZ + 30. 
The Manning parameter £m = ^BQrod/eL 30] for these 
rods is between 3.9 and 7.8, and counterion condensation 
is expected which, however, does not affect the shape 
of the rigid-rod macroion. We see again that the re- 
versed mobility for the cylindrical macroion is enhanced 
by the addition of small amount of monovalent salt, simi- 
larly to the strongly charged spherical macroion in Fig0 
The mobility for a ro d — 0.08e/a 2 is still reversed at 
n s i ~ 0.008/a 3 . Thus, the rod-shaped macroion is more 
persistent to monovalent salt than the spherical macroion 
of the same radius and surface charge density. The mobil- 
ity of the macroion with a ro d = 0.04e/a 2 is not reversed, 
whose surface charge density is twice large compared to 
that of the DNA. However, the polyelectrolyte counteri- 
ons can help to overcharge the macroion [3l], |32| • 

We compare in Fig^Jthe radial distribution functions 
of the spherical and cylindrical macroions of Fig(Sj and 
Fig|Sl under the same condition. The surface charge den- 
sity for these macroions is a ~ 0.08e/a 2 , radius R = 5a 
and salt ionic strength n s j ~ 0.006/a 3 . The peaks of the 
counterion charge density are of the same height for these 
cases, but the marked difference is that the monovalent 
counterions and coions are more closely concentrated to 
the cylindrical macroion than to the spherical macroion. 
Consequently, for the latter with non-reversed mobility, 
two dips below the neutrality occur in the trailing tail of 
the integrated charge profile in FigUJa). For the former 
with reversed mobility, there is only one and well-defined 
peak in the integrated charge profile of FigHJb), which 
reveals strong binding of ions to the cylindrical macroion. 
This case is associated with a well-developed lateral net- 
work of Z-ions on the surface of the macroion, as shown 
in the bird's-eye view Fig|5| Majority (70%) of the Z- 
ions are adsorbed on this rod macroion surface, which 
is more than sufficient for charge neutralization. Also a 
small fraction (20%) of the monovalent counterions are 
adsorbed to the macroion, unlike the case of the spheri- 
cal macroion of the same condition. Thus, the network 
is better packed by the counterions, which makes the 
surface ions less susceptible to thermal fluctuations and 
desorption. This supports more persistence of the rod- 
shaped macroion to monovalent salt than the spherical 
macroion of the same radius and surface charge density. 



V. SUMMARY 

In this paper, the charge inversion phenomenon was 
studied by electrophoresis, with the focuses on the effect 
of the asymmetric salt in size and valence, and the thresh- 
old of surface charge density of the macroion for both the 
spherical and cylindrical macroions. The criteria for the 
charge inversion obtained by the electrophoretic mobil- 
ity and the radial distribution functions of ions coincided 
except around the charge inversion threshold. 
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First, large coions compared to the counterions played 
a positive role in enhancing the charge inversion, while 
the coions with larger valence than the counterions did 
a negative role. More specifically, the reversed mobil- 
ity started at nearly null for very small coion radius, in- 
creased linearly with the ratio of the coion to counterion 
radii up to a~ /a + rs 1.5. Beyond this value, the mobil- 
ity decreased sharply due to large volume fraction of the 
coions and destruction of the complex by collisions. On 
the other hand, the mobility was reversed and largest for 
the monovalent coions, i.e. for the largest asymmetry of 
valences. The reversed mobility decreased with the ratio 
of the coion to counterion valences Z~ /Z + , and flipped 
to normal (non-reversed) for the coion valence exceeding 
that of the counterions. The mobility was reversed for 
divalent counterions and coions. 

Secondly, the addition of monovalent salt enhanced re- 
versed mobility for a strongly charged macroion at small 
ionic strength. This was due to global minimization of 
the electrolyte energy, and not only due to the local ad- 
justment of ions on the macroion surface. On the other 
hand, a large amount of monovalent salt screened the 
electrostatic interactions and suppressed charge inver- 
sion. For a weakly charged macroion, the enhancement 
regime was not detected and reversed mobility decreased 



monotonically with the increase in the monovalent salt 
concentration. 

Thirdly, the threshold of surface charge density under 
the monovalent salt was examined. The Z-ion (multiva- 
lent counterion) correlation energy at the threshold was 
several times of the thermal energy, Eq.J3J. A cylindrical 
macroion attracted counterions and coions more closely 
to its vicinity, formed a well packed network of the coun- 
terions, and was more persistent to the monovalent salt 
than the spherical macroion of the same radius and sur- 
face charge density. 
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